Collective excitations in unconventional charge-density wave systems 
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The excitation spectrum of the t-J model is studied on a square lattice in the large N limit in 
a doping range where a d-density-wave (DDW) forms below a transition temperature T*. Char- 
acteristic features of the DDW ground state are circulating currents which fluctuate above and 
condense into a staggered flux state below T* and density fluctuations where the electron and the 
hole are localized at different sites. General expressions for the density response are given both above 
and below T* and applied to Raman, X-ray, and neutron scattering. Numerical results show that 
the density response is mainly collective in nature consisting of broad, dispersive structures which 
transform into well-defined peaks mainly at small momentum transfers. One way to detect these 
excitations is by inelastic neutron scattering at small momentum transfers where the cross section 
(typically a few per cents of that for spin scattering) is substantially enhanced, exhibits a strong 
dependence on the direction of the transferred momentum and a well-pronounced peak somewhat 
below twice the DDW gap. Scattering from the DDW-induced Bragg peak is found to be weaker 
by two orders of magnitude compared with the momentum-integrated inelastic part. 

PACS numbers: 74.72.-h,74.50.+r,71.10.Hf 
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I. INTRODUCTION 

The electron density operator p on a lattice has in mo- 
mentum space the general form, 



= ^Y, T( k ' q)4(k + q)c CT (k). 

C kcr 



(1) 



c* a (k) , c<j (k) are electron creation and annihilation oper- 
ators with momentum k and spin projection a, N c is the 
number of primitive cells, q the transferred momentum, 
and k characterizes the relative motion of electron and 
hole. If both particles reside always on the same lattice 
site 7 is independent on k and p describes a usual den- 
sity wave. This wave may condense for some wave vector 
q different from the reciprocal lattice vectors of the un- 
modulated lattice and a conventional charge density wave 
(CDW) state is obtained. If the electron and hole explore 
different sites 7 depends on k and the condensation of this 
wave yields a CDW state with an internal symmetry de- 
scribed by the k dependence of 7. We will call this state 
in the following an unconventional CDW statei. Pro- 
posals for systems with an unconventional CDW state 
include organic conductors^ and the pseudogap phase 
of high-temperature superconductors2diSi&. 

Static and fluctuating density waves of the general 
form of Eq. Q can be observed by various experimental 
probes. The vector potential generated by the magnetic 
moment of neutrons influences the hopping elements of 
electrons via the Peierls substitution. Thus neutrons 
probe orbital magnetic fluctuations associated with den- 
sity fluctuations. As shown in Ref 3 the resulting neutron 
cross section is related to the imaginary part of a retarded 



density-density correlation function, where 7 is given by, 
Sine/ioteff 



7(k, q) 



fic|q| 



( — ( (1 — cos qj) cos kj + sin qj sin kj 
\ qj 



with 



ii x JL 
Mo \q\ 



(2) 



(3) 



The expression Eq. J5J applies to a solid of square layers 
with lattice constant a which we put, together with h to 
one in the following. t e f f is an effective nearest neighbor 
hopping integral, c the velocity of light, \x the magnetic 

moment of the neutrons, and po = J We neglect 
hopping between the layers as well as a small contribu- 
tion to 7 from second-nearest neighbor hoppings ~ t' . 
Similarly, inelastic light and non-resonant X-ray scatter- 
ing are determined by the Raman or the X-ray scattering 
amplitude which corresponds to the following expression 
for 7^, 



7(k, q) = 



a/3 



<9 2 e(k + q/2) 

' dk a dk/3 



(4) 



e l and e s are the polarization vectors of the incident and 
scattered light, respectively and e(k) the one-particle en- 
ergy. In contrast to RefsyQwe assumed that the Peierls 
substitution is made in the renormalized Hamiltonian 
so that the renormalized hopping t e ff and one-electron 
band appear in Eqs.© and 10} . For Raman scattering 
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the momentum transfer is practically zero whereas no 
such restriction exists for X-ray scattering. 

Microscopic models for interacting electrons usually 
contain interactions between local charge densities or 
spin densities, such as, the Coulomb or Heisenberg in- 
teraction, k-dependent densities of the form of Eq.lQJ 
become important if the exchange terms are competing 
or dominating the direct, Hartree-like terms. The exci- 
tonic insulator is such a case where the Coulombic ex- 
change terms not only create excitons but force them to 
condense into a new ground stated. These ideas were ap- 
plied to high-T c cuprates by Efeto\s^ assuming that the 
barely screened Coulomb interaction between layers pro- 
duces interlayer excitons and their condensation. More 
recently, it was recognized that the t-J model in the large 
iV limit (N is the number of spin components) shows a 
transition to a charge-density wave state with internal 
d-wave symmetry (DDW) at low doping^ with a tran- 
sition temperature T* . This model describes electrons 
in the CuOi planes of high-T c cuprates and identifies 
the pseudogap phase at low dopings with a DDW state. 
One may expect that near the phase boundary to the 
DDW state or in the DDW state density fluctuations of 
the form of Eq.QJ become important and, according to 
the above discussion, could be detected by neutron or X- 
ray scattering. It is therefore the purpose of this paper 
to calculate the general density response of this model, 
both above and below T*, and to make predictions for 
the magnitude and the momentum and polarization de- 
pendence of neutron and X-ray cross sections. 

The order parameter of a commensurate DDW state 
with wave vector (tt, it) is purely imaginary which follows 
from the hermiticity of the large N Hamiltonian. As a 
consequence, circulating, local currents accompagny in 
general density fluctuations in the DDW state produc- 
ing orbital magnetic moments localized on the plaquettes 
of the square lattice. The circulating currents fluctuate 
above and freeze into a staggered flux phase below T* . 
Most previous calculations considered q = quantities 
such as the frequency-dependent conductivityi2*ii*i2ii^ or 
Raman scatteringifi*A4. Calculations for finite momentum 
transfers^ to be presented below are interesting because 
they probe the local properties of the density fluctuations 
and the flux lattice. 

In section II we will first reformulate the large N limit 
of the t-J model in terms of an effective Hamiltonian 
which contains usual electron creation and annihilation 
operators (i.e., no constraint operators), renormalized 
bands and an interaction between several charge density 
waves. This effective model contains as a special case the 
models used in Refs lllllOl in discussing anomalous charge 
density waves. We then extend previous normal-state 
calculations of the density response, given by a 6x6 sus- 
ceptibility matrix \, to the DDW state using the Nambu 
formalism. The obtained expressions hold for all mo- 
menta q and are general enough to discuss the dispersion 
of collective modes and the momentum and polarization 
dependencies of experimental cross sections. 



Section III contains numerical results for the formulas 
presented in section II. In subsection A the density fluc- 
tuation spectrum in the d-wave channel will be discussed 
as a function of momentum, both above and below T* . In 
subsections B and C results for the cross sections for non- 
resonant inelastic X-ray scattering and for polarized and 
unpolarized neutron scattering will be given, and subsec- 
tion D compares the magnitude of q integrated neutron 
cross sections for spin and for orbital scattering. Section 
IV contains our conclusions. 



II. DENSITY RESPONSE IN THE D-CDW 
STATE 

In Refs lTr3ll7l it has been shown that the density re- 
sponse in the normal state of the t-J model at large N 
can be obtained by using the following two sets of density 
operators, 

pUq) = ^^S«(k,q)4(k+q)c CT (k), (5) 



^) = ]v-E^(k)4(k + qMk), (6) 

c ko- 

with 

£ Q (k,q) = (l,t(k + q) + J(q), 73 (k),74(k),7 5 (k), 76 (k)), 

(7) 

and 

FpQn) = (t(k),l,2J 7 3(k),2J 74 (k),2J 75 (k),2J 76 (k)). 

(8) 

Compared to Refs llrJT^ a slight change in the represen- 
tation of the basis functions has been made by using the 
symmetrized functions 

73.5 = (cosfea.±cosA; v )/2, 

7 4. 6 = (smk x ±s'mk y )/2, (9) 

where the subscripts 3 and 4 refer to the + and 5 and 6 to 
the — sign. Here, and in the following, we put the lattice 
constant of the square lattice to 1. i(k) and J(k) are the 
Fourier transforms of the hopping amplitudes ty and the 
Heisenberg term Jy , respectively. Including nearest and 
second nearest neighbor hoppings t and t', respectively, 
we have 

t(k) = 2i(cos k x + cos k y ) + At' cos k x cos k y . (10) 

The corresponding density-density Matsubara Green's 
function matrix is defined by 

Xa/j(q,»w») = - / ' dre-^ T < T T p^ (qr^qO) > . 
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The general solution for \ i n the normal state has been 
given in the large N limit and the dispersion of macro- 
scopic density fluctuations, determined by the element 
X12, have been discusse d 1 -948 . At lower temperatures and 
dopings the density component ^(k, Q) freezes in for a 
wave vector Q approximately equal to (w, tt), leading to 
the DDW state. 

The Hilbert space of the t-J model does not contain 
double occupancies of sites, its operators are therefore 
X and not the usual creation and annihilation operators 
c icr and Ci„ of second quantization. In the large iV limit 
the t-J model becomes, however, equivalent to the fol- 
lowing effective Hamiltonian in terms of usual creation 
and annihilation operators, 



(12) 



N 

H eff = £e(k)4 (k)c(k) - £ Xy a (q)pt(q 

ker Q = l q 

e(k) are the one-particle energies in the large N limit 

e(k) = 5 -t(k) - ^ • ± J2 cos(p x )/(e(p) - /i). (13) 
c p 

fi is the renormalized chemical potential and / the Fermi 
function. The second term in Eq. i|12fl represents an ef- 
fective interaction originating from the functional deriva- 
tive of the self-energy with respect to the Green's func- 
tion. This interaction is hermitean though its present 
form does not show this property explicitly because of the 
chosen compact form to represent it. Using Eqs.JSJ and 
© one easily verifies that the sum of the terms a = 3... 6 
represent the Heisenberg interaction, written as a charge- 
charge interaction and as a sum of four separable ker- 
nels. The terms a = 1,2 originate from the constraint 
which acts in H e ff as a two-particle interaction. H e ff 
and the original t-J Hamiltonian in the large N limit 
become equivalent if only bubble diagrams generated by 
the second term in H e ff and no self-energy corrections 
are taken into account. In the DDW state < p a (Q) > 
and < p' a (Q) > become non-zero for a — 5. We rewrite 
H ef f as, 

H,ff = ^ e (k)ct(k)c (T (k)-7V c <p^(Q)>pt(Q ) 



(14) 



a — 1 q 



The tilde at the density operators means that only their 
fluctuating parts should be considered. 

In order to be able to use usual diagrammatic rules in 
the DDW state we introduce the Nambu notation, i.e., 
the row vector iptO*-) — (c£ (k), c£(k + Q)) and the corre- 
sponding column vector obtained by hermitean conjuga- 
tion. We also make the convention that the momentum 
in ip or ip 1 are always taken modulo the reduced Bril- 
louin zone (RBZ), i.e., lie in the RBZ. The momentum 
sums over the (large) Brillouin zone (BZ) are then split 



into parts inside and outside of the RBZ and the two 
parts then combined using the Nambu vectors and Pauli 
matrices. We obtain then, 



P'M) = V4(k + q)4*(k,q)Vv(k), (15) 

c ktx 



Pf^ = ^E'^(k + q)i^(k)^(k), (16) 



with 



A(k) - i«(k)P (k + q,k)+i( 2 )(k).P e (k + q,k), 
F 2 (k) = P e (k + q,k), 

^(q) = ^(k)F (k + q,k), (17) 

for (3 = 3. ..6, 

Bi(k,q) = F 2 (k), 

Ba(k,q) = t (1) (k + q)P (k + q,k) 

+ (t( 2 >(k + q) + J(q))P e (k + q,k), 
B a (k,q) = F„(k)/2J, (18) 

for a — 3. ..6, and 

F (k + q,k) = p(k + q,k)CT 3 -i(l-p(k + q,k))a 2 , 
P e (k + q,k) = p(k + q,k)<ro + (1 - p(k + q, k))<7i. 

(19) 

tW(k) and ^ 2 )(k) are the Fourier transforms of the near- 
est and second-nearest neighbor hopping amplitudes, Co 
is the 2x2 unit matrix, and a\ , er 2 , 03 are Pauli matri- 
ces. The dash at the summation sign in Ea. l|l()l) means 
a restriction of the sum to the RBZ. p(k + q, k) with k 
in the RBZ is equal to one if k + q, reduced to the BZ, 
lies in the RBZ and is zero otherwise. Finally, the first 
and second terms in H e ff can be combined and read in 
Nambu notation as 

r(0) 



-e_(k)ff 3 + *75(k)cr 2 Vv(k). 



(20) 



The energies e± are defined by (e(k) ± e(k + Q))/2. <f> is 
an abbreviation for the amplitude of the DDW and must 
be real because of the hermiticity of H e f f . 

Using the above Nambu representation the suscepti- 
bility matrix Xa/3 is obtained by a bubble summation, 

xMl) =E^fe)(l+X (0, (9)) 1 | (21) 



Here we used the abbreviation q = (ii/ n ,q) with the 
bosonic Matsubara frequencies v n — 2nTn, where T is 
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the temperature. X (?) stands for a single bubble and 
is given analytically by the expression 



(?) = |-E' rr [ G ( fc +9)^( k 'q) G ( fc )^( k )]' (22) 

kn 



with the Nambu Green's function 

(iw„ + p — e+(k))er + e-(k)a 3 + $7 5 (k)f7 2 



G(k) 



(tw„-£7i(k))(iw„-^(k)) 



and the two eigenenergies 



(23) 



£ li2 (k) = e+(k) -/i± ^/e 2 _(k) +7l(k)<i> 2 . (24) 

After performing the trace over Pauli matrices and the 
sum over fermionic Matsubara frequencies in Eq. H22|) 
only the sum over momenta is left for a numerical evalu- 
ation. In the normal state we obtain 

(0)/ • \ 

— > E a (k, q)F p (k -— — — ; — . 25 

N c e(k + q) - e(k) - iuj n 

Two special cases of Ea. (|21|l should be mentioned. If 
q is parallel to the diagonal, i.e., of the form q = (q, q), 
the components 5 and 6 in the matrices, which are odd 
under reflection at the diagonal, decouple from the other 
components. At the points q = (0, 0) and (tt, it) they 
even decouple from each other so that Eq.JSTJ can be 
solved by division, for instance^, 



(0)/ 



(26) 



for q = (0,0) and q = {n,w). In the first case %55 
describes the B\ g component of Raman scattering^, 
in the second case amplitude fluctuations of the order 
parameter—. The other special case refers to a DDW 
without constrain1*i2ii£. Since the components 1 and 2 
were caused by the constraint they can be dropped in 
this case and Eq.J2U reduces to a 4x4 matrix equation. 
The four components arise because a nearest neighbor 
charge-charge interaction leads to four separable kernels. 
If this interaction is furthermore approximated by one 
kernel, associated with the order parameter, Eq. l|21|l re- 
duces to one scalar equation, and one arrives at the mod- 
els of Refs lHllol 

The above formulas allow to calculate density correla- 
tion functions which can be expressed by p' or p defined 
in Eqs.©-©. The most general case may, however, in- 
volve two additional densities po and p' in Eqs.JSJ and 
®), respectively, where Eq and Fq are general functions 
of k. Going over to the Nambu representation Eq and 
Fq are then linear combinations of the Pauli matrices 
and general functions of k. Using the diagramatic rules 
and the effective Hamiltonian one finds for the general 
susceptibility xoo(?), 



Xoo(?) 



XooH?) 



ExS2(9)Xa/j(9)xS } (9)- (27) 



5=0.077 
q=(7t,7t) 
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FIG. 1: Negative imaginary part of X55 of the t-J model at 
large N for t' = —0.3, r\ = 0.001, q = (tt,tv), and doping 
S — 0.077. Inset: DDW order parameter $ as a function of 
temperature T. 



Xoo (?) is given by Ea. 1)22(1 where both a and (3 have 
been replaced by 0. Similarly, the expressions for \^ 

and xfo are obtained from Ea. l22|) by replacing only a 
or (3, respectively, by 0. Xaf3 is given by Eci. (|2TJl . 



III. NUMERICAL RESULTS 
A. D-wave susceptibility X55 

From now on we will write t = —\t\ and use \t\ as the 
energy unit. Furthermore, the Heisenberg constant J and 
the doping S will be fixed to the values 0.3 and 0.077, re- 
spectively. The mean-field value for the superconducting 
transition temperature is then practically zero so that we 
deal with a pure DDW system. FigQ shows the negative 
imaginary part of the retarded susceptibility X55(q> w )> 
denoted by X55 m t ne following, as a function of lo in 
the normal state, using t' — —0.3 and a small imagi- 
nary part r\ = 0.001. We have chosen the component 
X55 and the wave vector q = (jr, tt) because the transi- 
tion to the DDW state occurs in this symmetry channel 
and near this wave vector. Decreasing the temperature 
from T = 0.032 to T = 0.028 shifts spectral weight in 
the overdamped spectrum from large to small frequen- 
cies until a sharp and intense peak appears very near 
to the transition point T* ~ 0.028. Decreasing further 
the temperature the spectrum becomes again broad and 
moves to higher frequencies. The inset of Fig^shows the 
temperature dependence of the DDW order parameter $ 
which exhibits the usual mean-field behavior. 

Fig[21 shows a three-dimensional plot of x'5'5 as a func- 
tion of u in the normal state at T = 0.039. The mo- 
mentum runs from the point T = (0, 0) to the points 
S = (tt, 7t) and X — (0, tt), and back to the point T. The 
curves for T, S and X are drawn as thick lines to facilitate 
the visiualization. According to Ea. (|25|) . and similar as 



FIG. 2: X5S m the normal state at T = 0.039 calculated 
for t' = -0.3, 8 = 0.077, and along the line T-S-X-T in the 
Brillouin zone. 

in the case of the usual charge susceptibility, X55 exhibits 
a singular behavior at small frequencies and wave vectors 
in the limit rj — > 0: Fixing the momentum to q = (0, 0) 
X5s(0,w) is zero for any finite frequency to. Putting to 
to zero and taking the limit q — * X55 approaches a 
finite value. This explains the absence of a thick line 
at r in FigEl Moving the momentum from r to S the 
zero-frequency peak of measure zero at Y aquires a finite 
width, i.e., a finite spectral weight, and at the same time 
peaks at a finite frequency. Both the width and the peak 
increase first rapidly with momentum. The peak posi- 
tion passes then through a maximum at u> ~ 0.04, starts 
to decrease, and reaches a minimum at S signalling the 
transition to the DDW state at lower temperatures. The 
width of the peak as well as the extension and intensity 
of the slowly decaying structureless background increases 
up to the point S. Moving along the direction S-X-T the 
low- frequency peak looses rapidly spectral weight, van- 
ishes practically at X, and recovers somewhat towards 
T. At the same time a well-pronounced, strong high- 
frequency peak develops, its frequency decreases mono- 
tonically along the above path, and its intensity shows 
a maximum near the point X. The dispersion of this 
high-frequency peak is ~ svak x between Y and X, i.e., 
it is identical with the sound wave due to usual density 
fluctuations^. This means that d-wave density fluctua- 
tions, described by X55, strongly couple to usual density 
fluctuations along the path S-X-T whereas such a cou- 
pling is forbidden by symmetry along T-S. 

FigEH shows the same 3d plot as Fig|2] but the calcu- 
lation is now performed at T=0 in the DDW state. Dif- 
ferences in the two figures thus have to be ascribed to 
the formation of the DDW. Taking the different scales in 
the two figures into account one sees that the intensity of 
the high-frequency peak due to usual density fluctuations 
did not change much. Big changes, however, occur at low 
frequencies. x'55 n0 longer vanishes at finite frequencies 
at the r-point but rather shows there a well-pronounced 
and strong peak. Moving away from the Y point this 
peak looses rapidly intensity whereas its frequency does 



FIG. 3: X55 at T=0. calculated for t' = -0.3, S = 0.077, and 
along the lines T-S-X-T in the Brillouin zone. 



not change much. 

To understand the origin and the properties of the low- 
frequency peak in more detail we have plotted in Fig0] 

X55 and X55 as a function of frequency for several q 
values between the Y and S point. The dashed line in 

(0)" 

the upper diagram represents X55 a ^ the wave vector 
q = (tt, 7t). It contains only interband transitions across 
the DDW gap, mainly between the points X and Y and 
between the hot spots on the boundaries of the RBZ. The 
intraband contribution is zero by symmetry at S. Since 
$ is about 0.07 the interband transitions lead to a broad 
peak somewhat higher than twice the DDW gap. The 
solid curve in the upper diagram of Fig0] represents x'55 
at S. It shows a much more pronounced peak than the 
dashed curve which originates from the denominator in 
Ea. (l26|) . i.e., which represents a collective excitation. It 
corresponds to the amplitude mode of the DDW, where 
the order parameter $ is modulated without changing 
its d-wave symmetry. Its frequency is somewhat smaller 
than 2<I> which is expected for a d-wave state. For an 
isotropic s-wave ground state the two frequencies would 
be exactly the same. 

The lowest panel in FigQ] shows the same susceptibil- 

(0)" 

ities at the Y point. X55 again contains no intraband 
but only interband contributions consisting of vertical 
transitions which probe directly the DDW gap. As a 

(0)" 

result, X55 exhibits a well-defined peak at 2$ and a 
long-tail towards smaller frequencies due to the momen- 
tum dependence of the order parameter. The solid curve, 
representing X55, exhibits a well-pronounced peak much 
below 2$. It corresponds to a bound state inside the d- 
wave gap created by multiple scattering of the excitations 
across the gap due to the Heisenberg interaction. Due to 
its large binding energy the density of states inside the 
gap is rather small at the peak position so that the width 
of the peak is substantialy smaller than that of the am- 
plitude mode at the momentum (n,ir). Previously^, we 
have associated this peak with the B\ g peak seen in Ra- 
man scattering in high-T c superconductors. Note that 
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FIG. 4: x's5 (solid lines) and Xss (dashed lines) for momenta 
along the diagonal q = (q, q) . 



FIG. 5: D-wave cross section for nonresonant, inelastic X-ray 
scattering for the t-J model at large N in the DDW state at 
T — 0. The momentum varies along the points T-S-X-T in 
the Brillouin zone. 



there is roughly a factor 5 difference in the scales of the 
top and bottom panels in FigQ] It means that the am- 
plitude mode of the DDW order parameter has a much 
larger spectral weight and also smaller width at the T 
than at the S point. 

For momenta q between T and S both intra- and in- 
to)" 

tcrband transitions contribute to X55 • The intraband 
part, reflecting the residual Fermi arcs in the DDW state, 
turns out to be in general substantially smaller than the 
interband part and consists of a broad peak which dis- 
perses roughly as sin(g) and which looses its intensity 
when approaching T or S. So only the low-frequency 
tails and the structure around the frequency 0.04 in the 
second panel from above can be attributed to intraband 
scattering, the remaining spectral weight is due to inter- 
band scattering. Going from r to S the solid lines show 
the dispersion of the amplitude mode which is in general 

fol" 

well below the main spectral peak in X55 , representing a 
resonance mode with a frequency and width determined 
mainly by the denominator in Eq. i|26[l . From the different 
scales used in each panel it is evident that the intensity 
of the amplitude mode strongly decays away from the T 
point in agreement with FigO 



B. Non- resonant inelastic X-ray scattering 

Specifying the polarization vectors of the incident and 
scattered light and using Eas. ((l3|l and (fTUfl for the elec- 
tron dispersion the cross section for inelastic X-ray scat- 
tering is determined by x'x an( i Xx is obtained as a linear 
combination of the susceptibility matrix elements Xap- 
As an example, consider the case e l = (1,1) /\/2, e s = 
(1, — l)/y/2, which yields for q — the B\ g component 
of Raman scattering. Choosing the symmetry direction 



q = [q, q), we find, 

2t 2 

xx(q,«w„) = — J^(cOS 2 (|)x55(q,^n) 

+ sin2 (|)X66(q,^n) - sin(g)x 56 (q,iu; n )), (28) 

with the effective nearest- neighbor hopping element t e f / 
determined by the large N dispersion of Eq.|^3J). The 
second-nearest hopping element t' drops out in the con- 
sidered symmetry. In Fig[S] we have plotted xx with- 
out this prefactor along the points T-S-X-T generalizing 
Ea. (l28|) to an arbitrary k point in the Brillouin zone. One 
feature of x'x ' IS that the high-frequency side bands due 
to usual density fluctuations are practically absent. This 
may be explained by the fact that the various contribu- 
tions to x'x have both signs causing large cancellation 
effects. Another difference between Figs|31 and occurs 
around the point S where additional structure is seen in 
Xx due to the component x&6- 

C. Inelastic neutron scattering 

Using the symmetry-adapted functions Eq.Q the ex- 
pression for 7, Eq.J2J), can be written as, 

7(k,q) = ^f^E^(^(k) ) (29) 

with 

B3.5 = — (1 - cosq x ) ± — (1 - cosq y ), (30) 

Qx q v 

-B4 6 = —sinq x ± —sinq v , (31) 

q x ' q y 

where the subscripts 3,4 refer to the + and 5,6 to the 
— signs, respectively. The charge correlation function, 




FIG. 6: x'n describing neutron scattering as a function of ui 
along the line T-S-X-T in the Brillouin zone in the DDW state 
at T = using t' = -0.3, S = 0.077, q z = 0.5, and rj = 0.05. 

describing neutron scattering, is given by, 

Xjv(q,^) = ^7 1 rzi J 2^ B ' 3 ( c l) jB5 ^)W' 5 ( t l' w )- 

ZJ C|q| 0,5 

(32) 

Fig. El shows Xjv(q, w) as a function of uj for momenta 
along the symmetry line T-S-X-r in the Brillouin zone for 
the fixed momentuma transfer q z = 0.5. The spin of the 
neutrons was assumed to be polarized along the z-axis. 
The spectra show essentially one more or less well-defined 
peak which disperses roughly between the frequencies 
0.1 — 0.2. Similarly as in Fig0] this peak describes a 
mainly collective excitation of the system related to a 
bound state inside the DDW gap. In spite of the rather 
large momentum transfer q z = 0.5 the peak intensity in- 
creases substantially towards the point T. This increase 
becomes more and more pronounced if q z is further low- 
ered until it diverges like l/|q| 2 in the limit q z — > 0. 

An important special case of Eq. l|32|) is scattering by 
unpolarized neutrons. Averaging over all spin directions 
for the neutrons yields the charge correlation function 
Xn,u for scattering with unpolarized neutrons, 

X2V,«(q, W) = ^ - F a / 9(q)Xa/3(q, w), (33) 

a/3 

with 

Poo = ^( 8 ^°| e// ) 2 < >«. ( 34 ) 

where < ... > av denotes an average over all directions of 
the moment of the neutron. We obtain, 

< S a (q)B /3 (q) > av = T^[{ql + q 2 z )f a ft3 

|q| 

~qxq y {f a gp + g a fp) + (ql + q 2 z )g a gp}, (35) 

with / 3 = / 5 = (1 - cos q x )/q x , g 3 = -g 5 = (1 - 
cosq y )/q y , f 4 = fe = sinq x /q x , g 4 = -g 6 = smq y /q y . 
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FIG. 7: x'n u describing scattering with unpolarized neutrons 
as a function of ui for momenta around the point F in the 
DDW state at T = using t' = -0.3, 6 = 0.077 and q z = 0.5. 

Xn,u(<\,u) diverges at small momentum transfers in the 
DDW state. To extract its singular part one may take 
the limit q — > in the /, g, and x functions. One obtains 
then, 

1 q 2 

|q| |q] 

In the normal state Xee vanishes at q = for any fi- 
nite frequency as can be seen from the explicit expression 
Ea. (|25|) . In contrast to that it is finite at q = in the 
DDW state describing a continuum of particle-hole ex- 
citations with zero total momentum across the gap. As 
a result xn,u diverges quadratically in the momentum 
qii parallel to the layers for q z = 0. This divergence is 
caused by the long-range part of the interaction between 
the magnetic moment of the neutron and the electrons 
in the layers. It is thus a real effect causing a singular 
forward scattering contribution in the cross section. Inte- 
grating x'n u over q diverges logarithmically in a strictly 
2d description at q z = whereas it remains finite in three 
dimensions. The large enhancement of x'n u near t ne 
point r is illustrated in Figd The curve at T shows a 
peak at ui ~ 0.1 well inside the gap of 2Ao ~ 0.14, which 
disperses to higher frequencies moving away from T, es- 
pecially, along the diagonal. More dramatic is, however, 
the rapid decay of the intensity of the peak away from 
r which can be seen from the different scales used in 
plotting the various panels. 

Another special case of Eq. (|32|l is obtained for an ar- 
bitrary polarization for the neutron moment and a van- 
ishing momentum transfer in the z direction, i.e., for 
q = (fej %> 0). The functions Bi assume in this case 
the form for i = 3, 4, 

Bi = —r-r(-q y fi + q x 9i), 

Molql 

B,+ 2 = —r-\{-qyfi - Qx9i)- (37) 

Moiqi 
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The polarization dependence of xn is simply cos 2 9 where 
9 is the angle between the neutron spin and the normal to 
the layers. xn assumes its largest value if the spin is per- 
pendicular and is zero if the spin is parallel to the layers. 
In this case xn diverges as l/|q| 2 for all polarizations of 
the neutron and the singular contribution is isotropic in 
the q z = plane. For q z ^= the dependencies of xn on 
the polarization and the momentum no longer factorize 
but can be worked out for the most general case using 
Eq.E2J. 



D. Comparison of scattering efficiencies from spin 
and orbital fluctuations 



As discussed in the previous sections the total magnetic 
neutron cross section contains, in addition to the domi- 
nating spin flip contribution, an orbital part. Above and 
not too far away from T* to the DDW state the scattering 
from fluctuating orbital moments produced by circulat- 
ing currents around a plaquette are mainly concentrated 
around the wave vector Q = (tt, it). Below T* there is an 
elastic Bragg component at wave vector Q and a fluctuat- 
ing part throughout the Brillouin zone. Though magnetic 
scattering from spin flip and orbital scattering have many 
features in common there are three main differences: a) 
In spin scattering the form factor reflects the electron 
density distribution within an atom whereas in orbital 
scattering the form factor is related to the plaquette ge- 
ometry which gives rise to the trigonometric functions in 
the functions Bp; b) The long-range dipole-dipole inter- 
action between neutrons and the electrons dominates in 
orbital scattering yielding a singular forward scattering 
peak below T c whereas the interaction between neutrons 
and electrons for spin scattering reduces to a local one; c) 
The intensity of orbital scattering is substantially smaller 
than that for spin scattering. For instance, a ratio of 1:70 
has been obtained in Refy for Bragg scattering. In the 
following we present results on frequency and in-plane 
momentum integrated scattering probabilities in the two 
cases. In this way the anomalous forward scattering con- 
tribution in orbital scattering can be taken into account 
in this comparison. We also recalculate the intensity for 
Bragg scattering and find it much smaller than the value 
calculated in Ref 3. 

The expression 2%^ u (q, + n(ui)) describes the 
scattering probability of an unpolarized neutron with 
momentum and frequency transfers q and u>, respec- 
tively, from orbital fluctuations. Integrating over fre- 
quency and the in-plane momentum qii the integrated 
scattering probability P or b becomes, using Ea. (|33l) and 
reinserting a and h for clarity, 



Por b ( q ,)=P B r r b a33 (q z ) + PZ n ( qz ) 



dyn. 



(38) 



2/3 



2/3 



/oo 
dw(l+n(w))x£ )t ,(q,w) 



qt|,a,/3 



In Eq. i|38|l P or b has been splitted into a static Bragg con- 
tribution P^rb" ano - a dynamic part P^b ■ Correspond- 
ingly, the tilde in Eq. Q40(l indicates that only the fluctuat- 
ing part of the densities should be used in the correlation 
functions. 

Similar considerations apply to unpolarized neutron 
scattering from spin fluctuations. The relevant suscep- 
tibility in the paramagnetic state ia22, 



-7-1 — Xz,(q,^)- 

D V mca* / 



(41) 



F(q) is the atomic structure factor and x zz (q, w) the zz 
component of the electronic spin susceptibility. Using the 
fluctuation-dissipation theorem, approximating F(q) by 
F(0) = 1, and assuming no magnetic coupling between 
the layers, the total scattering probability P s by spin 
fluctuations becomes 2 ^, 



12 \mca i J 



(42) 



Most interesting are the ratios of integrated orbital to 
spin scattering, i.e., R Bra aa = P Br b a " jP s and R dyn = 
Porb 1 /Ps- Using the above results one obtains, 



R dyn {q z ) 



A- 



N, 



— T 

2/3 L^i 



(B a (q)Bp(q)) c 
(aq) 2 



c q||,a>/3 

d(h)u> (l + nH) X ^(q/|f|,o;), (43) 



with the dimensionless constant, 
4 



tt(1 - S)(J/\t\) \H 2 /(ma 2 ) 



(h 2 /(ma 2 )) ' (44) 



and 



R B ^(q z ) = A { -^0^(p' 5 (Q))(plm- (45) 

Using \t\ = 0.5 eV we obtain from Ea. (fT3^ |t e //| ~ 
0.07 eV, which yields togther with a = 3.856 A A - 0.09 
for our doping S — 0.077. Comparison of Ea. (|14|) with 
Eq.(j23) gives 

(^(Q)X/4(Q)> = ~ 0.082, (46) 

for T = and S = 0.077. Furthermore, from Ea. (|35|) we 

get 



Por r b a "(<lz) = 2F 55 (Q,qM(Q))(pUQ)) ) (39) 



(47) 
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FIG. 8: Ratios R Bva aa an d R d v n for elastic and inelastic 
orbital to spin scattering, respectively. 

which yields 

i? B ™ff(0) = 3- 1(T 5 . (48) 

Unlike in the case of spin scattering the momentum and 
frequency sums in Ea. l|43|) are non-trivial and we have 
carried them out by numerical methods. 

Fig. 8 shows the results for R d v n and R Bra aa as a 
function of q z at zero temperature. As expected R d v n 
increases strongly at small q z due to its logarithmic sin- 
gularity at q z = 0. Integrating R d y n also over q z yields a 
value of about 10 -2 which is more than 2 orders of magni- 
tude larger than R Bra ^. The reason for this large differ- 
ence is due to the fact that in the momentum integration 
in R d y n mainly the small momenta contribute where 1 /q 2 
is larger than 1. In contrast to that, R Bra ss probes the 
integrand of R dyn at the large value |Q| where 1/q 2 is 
very small. The extreme small value R Bra ss ^ 10~ 5 casts 
doubts on the identification of the observed Bragg peak 
at Q with the Bragg peak due a DDW—. On the other 
hand, if a DDW state is realized in some material we 
think that the resulting strong forward scattering peak 
in neutron scattering would be observable. According to 
our predictions it would exhaust roughly 3 per cent of the 
sum rule for spin scattering and thus be of the same order 
of magnitude as the recently detected elastic scattering 
from weak magnetic moments directed at least partially 
perpendicular to the layers^ in high-T c superconductors. 

IV. CONCLUSIONS 

The excitation spectrum of the t-J model at large 
N was studied for dopings where a DDW forms below 



a transition temperature T*. The density response in- 
cludes "conventional" local density fluctuations, charac- 
terized by the energy scale t, where practically all spec- 
tral weight is concentrated in a collective sound wave 
with approximately a sinusoidal dispersion. In addition, 
and this was the main topic of this investigation, there 
is a "unconventional" density response from orbital fluc- 
tuations caused by fluctuating circulating currents above 
T* and fluctuations around a staggered flux phase below 
T*. The corresponding spectra are again mainly collec- 
tive, their energy scale is J or a fraction of it, and they 
can be viewed as order parameter modes of the DDW. At 
small momentum transfers the spectra are dominated by 
one well-defined peak corresponding to amplitude fluctu- 
ations of the DDW. For larger momentum transfers, un- 
conventional density fluctuations tend to become broad 
and their intensity suppressed. This is especially true 
near the point X where practically no weight is left for or- 
bital fluctuations at low frequencies, but instead a rather 
sharp peak due to sound waves appear at large frequen- 
cies reflecting the coupling of the two kinds of density 
fluctuations. 

In principle the above orbital fluctuations show up 
in the Raman, X-ray and neutron scattering spectra. 
The predicted scattering intensities are, however, rather 
small. The reason for this is that, on the one hand, a 
rather small effective nearest-neighbor hopping element 
t e ff is needed to provide a sufficiently large density of 
states and the instability to the DDW, on the other hand, 
the square of t e ff appears as a prefactor in the various 
cross sections due to the Peierls substitution. As a re- 
sult, the q integrated inelastic cross section is only a few 
per cents of that for spin scattering and the elastic scat- 
tering from the Bragg peak is even much weaker. From 
our calculations it is also evident that a DDW state in 
the unconstrained, weak-coupling case would have very 
similar properties as in the t-J model. In particular, the 
curves for R d y n and R Bra 39 would closely resemble those 
in Fig|S| An experimental verification of a DDW state 
via inelastic neutrons seems to be not unfeasible because 
a good deal of the scattering intensity is concentrated in 
a well-pronounced, inelastic peak somewhat smaller than 
the DDW gap and in small region around the origin in k 
space. 
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